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Abstract. Two dimensional Fourier spectra of near- 
infrared images of galaxies provide a powerful diagnostic 
tool for the detection of spiral arm modulation in stel- 
lar disks. Spiral arm modulation may be understood in 
terms of interference patterns of outgoing and incoming 
density wave packets or modes. The brightness along a 
spiral arm will be increased where two wave crests meet 
and constructively interfere, but will be decreased where 
a wave crest and a wave trough destructively interfere. 
Spiral arm modulation has hitherto only been detected 
in grand design spirals (such as Messier 81). Spiral arm 
amplitude variations have the potential to become a pow- 
erful constraint for the study of galactic dynamics. We il- 
lustrate our method in two galaxies: NGC 4062 and NGC 
5248. In both cases, we have detected trailing and lead- 
ing m=2 waves with similar pitch angles. This suggests 
that the amplification mechanism is the WASER type II. 
In this mechanism, the bulge region reflects (rather than 
refracts) incoming waves with no change of pitch angle, 
but only a change of their sense of winding. The ratio 
between the amplitudes of the leading and the trailing 
waves is about 0.5 in both cases, wherein the higher am- 
plitude is consistently assigned to the trailing (as opposed 
to leading) mode. The results are particularly significant 
because NGC 5248 is an optically grand design galaxy, 
whereas NGC 4062 is optically fiocculent. NGC 4062 rep- 
resents the very first detection of spiral arm modulation 
in the stellar disk of an optically flocculent galaxy. 

Key words: Galaxies: spiral - Galaxies: structure - 
Galaxies: kinematics and dynamics - Galaxies: individual 
(NGC 4062, NGC 5248) - Methods: numerical 



1. Introduction 

Within t he fra mework of the Density Wave Theory (Lin 
and Shu 1964 ), the spiral arms of disk galaxies are the 



manifestation of sets of travelling waves. In the presence of 
both leading and trailing sets of waves, the modal theory 
of galactic spiral structure (Berlin et al. 1989a, 1989b) pre- 
dicts that the amplitude of a grand design two-arm spiral 
pattern will oscillate with radial distance from the center 
because of the interference of wave packets or modes which 
are propagating inward and outward, being reflected off a 
central bulge (Lin 197C| ; Lau et al. 
1981. 
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The swing amplification mechanism (Toomre 1981) 
may be responsible for the appearance of the modes them- 
selves. The existence of trailing and leading waves is in- 
ferred by looking at the unmistakable interference pat- 
terns on the density contours along spiral arms (see, e.g., 
Toomre's Figures 10 and 12). The presence of interference 
patterns is betrayed in the discontinuity of the contours, 
or in other words, in the modulation of the density (or 
surface brightness) along each arm. The brightness along 
a spiral arm will be increased where two wave crests meet 
and constructively interfere, but will be decreased where 
a wave crest and a wave trough destructively interfere. 



In the swing amplification models of Toomre (1981), 
the absence of a bulge or of a central concentration 
makes the modelling unrealistic. The modes which grow in 
Toomre's models are fast evolving. In contrast, the modal 
theory of galactic spiral structure assumes the presence 
of a central bulge. These bulges act as reflectors (often 
termed the 'Q-barrier'), resulting in a quasi-stationary 
modal pattern (Thomasson et al. 199C; Elmegreen and 



Thomasson |1993t Fuchs |1991| , ^OOO] ). Clearly, such mod 
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els are much more appropriate to the dynamics of spiral 
galaxies. 

Several mechanisms have been proposed for the main- 
tenance of spiral structure in galaxies. The usual process of 
wave propagation, with feedback and over-reflection, can 
maintain the wave pattern. Mechanisms such as turbulent 
dissipation and shock formation in the gaseous Population 
I component can also play a role in the self-regulation of 
spiral modes (Berlin et al. 1989a, 1989b). 
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Symmetric spiral arm amplitude modulations indica- 
tive of underlying wave modes have hitherto only been 
detected in the grand design galaxies M51, M81, MlOO 
(Elmegreen et al. 1989D and in the multiple arm galaxy 
MlOl (Elmegreen 1995| ). Arm variations in two other 
galaxies were discussed by Grosb0l ( |1988D . M81 is one of 
the best studied spiral galaxies, and the amplitude data is 
very useful to constrain the model parameters within the 
modal theory. The m odal theory has been applied to M81 
by Lowe et al. ( 1994 ), wherein the observed arm modula- 
tion was modeled. 

Near-infrared images reveal the old stellar Population 
II disk component of spiral galaxies^, while optical images 
show the rich variety of responses of the young Population 
I component to the underlying older stellar population 
responsible for the dynamics of the galaxies (Frogel et al. 



Table 1. Parameters of the galaxies 



1996). The young Population I disk component may only 



constitute 5 percent of the dynamical mass of the disk of 
a galaxy. For studying mass distributions of disk galaxies, 
near-infrared images are essential (Block and Wainscoat 



1991; Block et al. 1994; Quill en et al. 1994, 199^; Frogel 



et al. |l996t Block and Puerari |1999[ Block et al. pOOq) . 

Optically thick dusty domains in galactic disks can 
completely camouflage or disguise underlying stellar struc- 
tures. Dust extinction is highly effective whether or not 
the dust lies in an actual screen or is well intermixed 



with the stars (Elmegreen and Block 1999). The pres- 
ence of dust and the morphology of a galaxy are inextri- 
cably intertwined: indeed, the morphology of a galaxy can 
completely change once the Population I disks of galax- 



ies are dust penetrated (e.g.. Block and Wainscoat 1991 



Block et al. 1994). Dust can completely obscure two armed 



grand design structure in some optically flocculent galax- 



1999) 



ies (Thornley 1996 Block and Puerari 199£ ; Elmegreen et 
al. 



In this paper, we propose a morphological method, 
based on the bi-dimensional Fourier transform, to detect 
the existence of structures with a different winding sense 
(trailing and leading patterns) in the same galaxy. The 
galaxies for which spiral arm modulations have hitherto 
been detected (e.g., M81) are nearby. The Fourier spectra 
offer an unambiguous way of identifying both leading and 
trailing wave packets in galaxies which are not restricted 
to be relatively close; the method can be applied to any 
spiral whose stellar spiral arms are resolved. 



^ Observational astronomers invariably restrict the terminol- 
ogy of 'Population IF for the stars in the halo of a galaxy, and 
refer to the 'young Population I disk' and the 'old Population I 
disk'. When modelling the disks of galaxies, however, theorists 
find it convenient to distinguish the two dynamically different 
gaseous and stellar components by 'gaseous Population I disk' 
and 'evolved stellar Population II disk', and we retain that 
terminology here. It must, however, be emphasized that by an 
'old stellar Population II disk' we are not referring to any true 
metal poor Population in the halo. 



NGC 4062 NGC 5248 



Type-^ 

DP class*" 
arm class*^ 

Filter 



Sc(s)II-III 
-19.44 
E/3 
3 
100° 
64° 
H {1.65fj,m' 



" Sandage and Tammann 1987 
Block and Puerari 



1999 



Sbc(s)I-II 
-21.19 
E/3 
12 
110° 
43° 
K' (2.1/im) 



'^ Elmegreen and Elme green 
de Vaucouleurs et al. 1991 



1987 



The Fourier method is applied to the near-infrared im- 
ages of two galaxies which optically could not be more 
different: one is flocculent (NGC 4062) whereas the other 
(NGC 5248) is grand design. The morphological appear- 
ances of NGC 4062 and NGC 5248 in the dust penetrated 
regime are carefully discussed below. 

2. Data and analysis 

The H band {1.65fim) image of NGC 4062 is part of the 
OSU (Ohio St ate U niversity) Bright Spiral Galaxy Survey 
(Frogel et al. |l996| ; Eskridge et al., in prep.). NGC 5248 
was observed in the infrared (K' 2.1fim) at the Observa- 
torio Astronomico Nacional at San Pedro Martir, Mexico, 
and forms part of a larger project on deep K' imaging. 
Details of reduction will be given elsewhere. 

Since our focus is morphology, no calibration frames 
are required here. Using the IRAF| task IMEDIT, the 
images were cleaned of any foreground stars. The galax- 
ies were then deprojected using the IRAF ROTATE and 
MAGNIFY routines. Deprojection parameters (position 
angle PA and inclination w) are listed in Table |^; also 
given in that Table are van den Bergh luminosity classes 
and blue absolu te ma gnitudes as determined by Sandage 
and Tammann ( |l987j) . Morphological parameters such as 
DP and arm classes are discussed below. 

Once the galaxies are corrected to a 'face-on' orienta- 
tion, we applied the program 2dfft (see the appendix of 
Schroeder et al. 1994), which calculates the fast Fourier 
transform of a given image using a basis of logarith- 
mic spirals. As shown in an extensive work by Danver 



(1942), logarithmic spirals appear to be the best math- 
ematical description for galactic arms. In a more recent 



work, Kennicutt (1981) concluded that logarithmic, as 
well as hyperbolic spirals, are good representations of 



^ The IRAF package is written and supported by the IRAF 
programming group at the National Optical Astronomy Ob- 
servatories (NOAO) in Tucson, Arizona. NOAO is operated 
by the Association of Universities for Research in Astronomy 
(AURA), Inc. under cooperative agreement with the National 
Science Foundation (NSF). 
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galactic spiral arms. As discussed elsewhere (e.g., Con- 
sidere and Athanassoula 1982; Puerari and Dottori 1992D , 
the choice of logarithmic spirals docs not constrain the 
analysis. They only form the basis in a vector space for 
the decomposition. If only a few coefficients - typically, 
one or two - are required to reconstruct the original im- 
age using inverse Fourier transforms (as is the case here), 
the choice of logarithmic spirals is indeed appropriate. 

The Fourier method has been extensively discussed in 
a number of paper s (e.g., Kaln ajs [1975 ; Considere and 
Athanassoula |l982|; lye et al. |l982[~ Puerari and Dot- 



ton 



1992 ; Puerari 1993 , amongst others). In the Fourier 
method, an image is decomposed into a basis of logarith- 
mic spirals of the form r—roexp{—^9). The Fourier coef- 
ficients A{p, m) can be written as 



A{p,m) — — / / I{u,0)exp[—i{m9 + pu)]dud9 

Here, u = In r, r and 9 are the polar coordinates, 
m represents the number of the arms, p is related to the 
pitch angle P of the spiral by P=atan(— m/p), and /(u, 9) 
is the distribution of light of a given deprojected galaxy, 
in a (In r, 9) plane. _D is a normalization factor written as 



D 



-\-OQ 



I{u,9)dud9 



In practice, the integrals in w = In r are calculated 
from a minimum radius (selected to exclude the bulge 
where there is no information of the arms) to a maximum 
radius (which extends to the outer limits of the arms in 
our images). 

The inverse Fourier transform can be written as 



Siu,9)=J2Srniu)e 



irn6 



where 



Smiu) 



D 



e^"47r^ 



G^{p)A{p,m)e'P" dp 



and G„i{p) is a high frequency filter used to smooth the 
A(p, to) spectra at the interval ends (see Puerari and Dot- 
1992D , and it has the form 



tori 



Gmip) = exp 



1 /„ _ „™ 

2 \ 25 



where p^ax is the value of p for which the amplitude of the 
Fourier coefficients for a given m is maximum. The chosen 
interval ends (p+=-|-50 and p_=— 50), as well as the step- 
size dp=0.25, are suitable for the analysis of galactic spiral 
arms. 

In Table ^ we present the values for the dominant to=2 
components in the Fourier spectra. In that Table, Pmax is 



Table 2. Values of the to=2 Fourier spectrum 



NGC 4062 NGC 5248 



Pmax 


-5.75 


-3.75 


Pmax 


4.5 


3.5 


\A{pLa.,'i)\ 

|A(p^„„2)| 
\A^\/\A^\ 


2.1E-3 


1.5E-1 


1.2E-3 


8.0E-2 


0.54 


0.54 




19° 


28° 


pL 


-24° 


-30° 


$^ 


-34° 


-87° 




-75° 


-72° 



the value where the spectrum for m=2 peaks, \ A(j>max, 2) 
is the amplitude of the peak, and denotes the 

corresponding ratio between the leading (L) and the trail- 
ing (T) amplitudes. P is the pitch angle of the spiral, re- 
lated to Pmax by P=atan(— 2/pmaa:)- ^ is the phase of the 
spiral arm, calculated as ^—a.tan{Im[A]/Re[A]), where 
/to [A] and i?e[A] are the imaginary and the real part of 
A{prnax, 2), respectively. A diagram explicitly showing the 
definition of trailing as opposed to leading spiral arms may 



be seen in Figure 15 of Athanassoula (1984), in Figure 2.8 
of Bertin and Lin (1996) or in Figure 6.5 of Binney and 
Tremaine ( |1987D . 

It is evident from Table || that the pitch angle for the 
trailing and the leading components for each galaxy is very 
similar. A classification scheme of spiral galaxies in the 
near-infrared was recently proposed by Block and Puerari 
( |1999D . Galaxies are binned into three groups a, /? and 7 
based on the pitch angle of the arms, robustly determined 
from Fourier spectra. Even-sided (as opposed to lopsided) 
galaxies have a dominant to=2 component in the Fourier 
spectra, and they are designated in this scheme by an 
'E'. NGC 4062 and NGC 5248 both belong to the dust 
penetrated E/3 class (Table |l|) . 

The ratio between the amplitudes of the leading and 
the trailing patterns is also the same for the two galax- 
ies. This is particularly interesting, since there is a large 
difference in the linear size of the two galaxies (adopting 
Ho=50 km/sec/Mpc, the linear diameters of NGC 4062 
and NGC 5248 are 18 and 42 kpc, respectively). Note fur- 
thermore that the ratio of amplitudes does not depend on 
the absolute magnitude of the parent spiral, neither on 
the arm class (see Table |^) . 

In the determination of Fourier coefficients and of pitch 
angle, careful deprojections to face-on are of course nec- 
essary. Mean uncertainties of position angle and inclina- 
tion angle as a function of inclination are drawn in Fig. 



2 of Considere and Athanassoula ( |1988| ). For NGC 4062, 
we find (using eight deprojected runs, wherein inclination 
and position angle are systematically varied) that incor- 
rect deprojections can introduce a maximum difference of 
13% in the reported ratio (see Table |). For 

the pitch angle, the maximum difference we find is only 
5°. The situation is slightly more complex for NGC 5248, 
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Fig. 1. Left: Deprojected H image of NGC 4062. The vertical bar represents 5 kpc (Ho=50 km/sec/Mpc). Right: The 
contours of the inverse Fourier transform for the m=2 mode are overlayed on the deprojected image. The contours are 
the real part of the complex spatial function S2{u, 9)=S2{u)e^'^^ . The lowest and the highest plotted contours represent 
3% and 30% of the maximum amplitude for the m=2 mode, respectively. 



because of its smaller inclination (and hence, larger uncer- 
tainties in the deprojection angles) and its asymmetry. For 
a few incorrect deprojections (mainly when we deproject 
the image with both an incorrect PA and an incorrect to), 
the leading component does not appear clearly. In other 
cases, the errors in the |A^|/|v4^| ratio and in the pitch 
angles are of the order of the errors listed for NGC 4062. 
The results are fully consistent with the findings of Block 
et al. (2000) wherein determination of pitch angles from 
Fourier spectra are found to be surprisingly robust; galax- 
ies do not move from one dust penetrated (DP) class to the 
next. NGC 4062 and NGC 5248 remain dust penetrated 
(3 class. 



2.1. NGC 4062 

Optical images of the galaxy NGC 4062 (Sandage and 
Bedkc 1994, plate 265) reveal numerous patches of star 
formation with no grand design spiral structure from a 
density wave in the underlying stellar disk. Such patchy 
structure occurs in over sixty percent of isolated, non- 
barred galaxies, giving them a flocculent, fleece-like ap- 
pearance (Elmegreen and Elmegreen 1987). A character- 



NGC 4062 belongs to arm class 3, described as 'frag- 
mented arms uniformly distributed around the galactic 
center' (see Elmegreen and Elmegreen 1987 ). What is 
so striking in the dust penetrated, near-infrared regime 
is that NGC 4062 presents a remarkable, bisymmetrical 
grand design morphology (see Figure |l]): two arms span- 
ning over 90 degrees in azimuth. It is evident that the 
young Population I and old stellar Population II disks of 
NGC 4062 actually decouple|. 

Another surprise is that the Fourier spectra of NGC 
4062 do not present a single peak for the m=2 mode (Fig- 
ure H) . The presence of two peaks betrays the existence of 
two different spirals. One would have the form of an 'Z' 
(coming from the peak with maximum at p < 0) and the 
other having a 'S' form. This situation is exactly what is 
expected within the framework of the swing amplification 
theory or in the modal theory of galactic spiral structure: 
a strong trailing pattern and a weaker leading one. 



istic of flocculent galaxies such as NGC 4062 is that the 
optical patches, by definition, span only a small range in 
azimuth. 



What is meant by decoupling is that a galaxy may show 
two different morphologies when examined optically and in the 
near-infrared regime. For example, NGC 309 is classified as Sc 
in the optical, but ap pears as a SBa in the near-infrared (Block 
and Wainscoat 1991). Two different morphologies in the same 
galaxy may co-exist via a f eedba ck mechanism or dynamical 
thermostat (Bertin and Lin f996| ). 
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Fig. 3. Left: Deprojected K' image of NGC 5248. The vertical bar represents 5 kpc (Ho=50 km/sec/Mpc). Right: The 
contours of the inverse Fourier transform for the m=2 mode are overlayed on the deprojected image. The lowest and 
the highest plotted contours represent 8% and 90% of the maximum amplitude for the bisymmetric mode, respectively. 



' 1 ' 1 ' 1 ' 1 
: N4062 m=l : 


'1 


m=2 7 








7 m=5 7 


■-i — r^r 1 


m=6 



-40 -20 20 40 -40 -20 20 40 
P P 

Fig. 2. Fourier spectra for NGC 4062. Note the absence of 
odd modes which is indicative of a highly bisymmetrical 
light distribution. 

2.2. NGC 5248 

Grand design spirals (where two dominant arms may span 
many degrees in azimuth) are completely different from 



flocculent galaxies. Possibly two of the best studied grand 
design sphals are Messier 81 (NGC 3031) and M51 (NGC 
5194), where two long symmetric arms dominate their 
optical disks. Elmegreen and Elmegreen (1987) devote 
classes 9 through 12 to the grand design bin. 

NGC 5248 is a magnificent, grand design optical spec- 
imen. Its arm class is 12 (the same class to which M81 
and M51 belong). In grand design galaxies, density waves 
are believed to have organized young stellar associations 
to fo rm the symmetric optical spiral pattern (Elmegreen 
1995) and one often finds a strong coupling between the 
young Population I and the older stellar Population II 
disks. This is the case for NGC 5248. 

This galaxy, however, is not quite as symmetric as 
NGC 4062 is in the near- infrared K' image (see Figure ||). 
The two trailing arms have a small difference in their wind- 
ing angle, with one arm being slightly more open than the 
other one. In the Fourier spectra, m—1 components reveal 
this asymmetry, and this is attested to by the relatively 
large m=l component of NGC 5248 (Figure 

Nevertheless, the Fourier spectra in Figure ^ reveal 
something quite remarkable: the old stellar Population II 
disk of an optically grand design galaxy, NGC 5248, can 
be almost identical to that of an optically flocculent, NGC 
4062 (compare Figures | and |). NGC 5248 also has two 
peaks for the m—2 mode, one at p < and another for 
p > 0, which proves the existence of two wave-trains with 
different winding senses: one trailing, the other leading. 
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Fig. 4. Fourier spectra for NGC 5248. The peak at m^l 
represents the asymmetry of the galaxy (see Figme ^. 
Note the striking similarity between the spectrum for m=2 
in this Figure and that in Figure ||. 

2.3. The trailing and the leading spirals 

In order to shed further information about the modes we 
have detected, we calculated the inverse Fourier transform 
in a different way to that illustrated in Figs. and I The 
goal is to separate the trailing and the leading compo- 
nents, and in order to study their individual characteris- 
tics, we calculate the 5*2 (u) functions by considering only 
p < OT p > 0. In other words, for the trailing mode, we 
use 



D 



G^ip)A{p,m)e'P" dp 



and for the leading one. 



D 



o2u 



47r^ 



p+ 



G„,{p)A{p,myP^ dp 



In Figs. H and ^ we plot the 5'|'(m) and the func- 
tions for NGC 4062 and NGC 5248, respectively These 
functions peak where we find the first maximum in the 
arm/interarm contrast (see Figs. ^ and |l^ below). Also 
plotted in Figs. |^ and ^ is the ratio between the ^^^(w) 
and the Sjiu) functions. This ratio is important for com- 
parisons with theoretical studies. As one can see, the ratio 
is of the order of 0.4—0.5 for almost all the spatial extent 
of the spiral arms in both galaxies. The ratio shows an 
increase towards smaller distances (nearer to the bulges 



N4062 




trailing . 
leading ^ 
ratio o 



< 



radius (kpc) 

Fig. 5. The surface density S'|'(r) and S'J(?') functions 
of NGC 4062, and the ratio S^/Sj. The cuts in radius 
are the minimum and maximum radii used in the Fourier 
analysis. 



of the galaxies) and also for larger radii (probably near to 
the CO- rotation radius - see below). 

By using the 5'|'(u) and the <S'J(u) functions, we can 
separate the leading and the trailing spirals of each galaxy. 
In Figs. and || we plot th e phase of each component (see 
Puerari and Dottori |l997|) for NGC 4062 and NGC 5248, 
respectively. The phase represents the azimuthal position 
of the maximum intensity of each spiral. So, when the dif- 
ference between the phases of the leading and the trailing 
components is 0° or 180°, the spiral pattern of the galaxy 
must show a maximum. In contrast, when the phase dif- 
ference is 90°, the spiral arms of the galaxy must show a 
minimum. For NGC 4062, we find 2.2 and 3.9 kpc for the 
positions of the maxima, and 2.9 kpc for the minimum. 
These values are in complete agreement with the peaks 
and the dip in the arm/interarm contrast (see Fig. ^ be- 
low). The situation is a little more complex for NGC 5248. 
For this galaxy, the values for the positions of the maxima 
are 3.6 and 8.3 kpc, and for the minimum, we find 5.2 kpc. 
These values do not fit the arm/interarm contrast as well 
as does the data for NGC 4062 (see Fig. below). This 
discrepancy could be explained by the fact that NGC 5248 
is not quite as symmetric as NGC 4062. 

Future models might be able to use these leading and 
trailing waves to locate the co-rotation radius. Leading 
modes propagate outward inside co-rotation and inward 
outside co-rotation, so that a leading spiral commenc- 
ing from the nuclear region would only extend to the 
co-rotation radius; in other words, co-rotation would be 
where the leading spiral ends (C. Yuan, private commu- 
nication). The difficulty in the present study in placing 
the co-rotation resonance is that the leading arm is not 
directly seen in our near-infrared images; rather, its exis- 
tence is inferred from the arm modulation. 
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Fig. 6. Same as Fig. |, but for NGC 5248. 
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Fig. 8. Same as Fig. 0, but for NGC 5248. 



2.4- Arm/Interarm contrast 

It is well known that the arm/interarm contrast in spi- 
ral galaxies usually increases with radius (Elmegreen and 
Elmegreen 1984 ) . While some of the galaxies in their sam- 
ple revealed arm/interarm profiles with a simple sinusoidal 
pattern in the I band (0.85^to), a large percentage of them 
showed a chaotic behaviour. Some galaxies can be opti- 
cally thick even at I. We have calculated the ar m/inte rarm 
contrast following Elmegreen and Elmegreen (1984). We 
have drawn azimuthal profiles for a number of radii, and 
by using an interactive plotting program, we have taken 
a mean value at the maxima (where the spiral arms are 
located) and a mean value at the 'troughs' (the interarm 
regions) . The final value is calculated from intensity values 
as follows: 



arm/interarm — (arm — interarm)/0.5(arm + interarm) 

In Figures ^ and |l^ we plot the calculated 
arm/interarm contrast for NGC 4062 and NGC 5248, 
respectively. The maximum values for NGC 4062 (the 
flocculent galaxy) are approximately 3 times less than 
those for NGC 5248 (grand design). As one can see, the 
arm/interarm contrast increases with radius, but not in a 
monotonic way. The modulation of the intensity is caused 
by the interference between the incoming and the outgoing 
density waves and is clearly evident in our plots. The sinu- 
soidal behaviour is even more remarkable in the flocculent 
NGC 4062. 

It is important to note that the values for the 
arm/interarm contrast calculated for these two galax- 
ies are in complete agreement with other studies (eg., 
Elmegreen et al. 199(31 ). Values of arm/interarm contrast 
as high as 2 are not unreasonable (see Fig. 6 of Elmegreen 
199^ A near-infrared study of MlOO by Gnedin 



et al. 

et al. (1996) yields an arm/interarm contrast of 3. High 
arm/interarm contrasts do not imply unrealistic high per- 
turbations of the velocity fields. Rather, what is important 
is the ratio of the arm mass to the enclosed galaxy mass at 
a given radius. The arm is but a small fraction of the en- 
tire mass of the galaxy inside any given radius, so that the 
streaming motions are small, even for high arm/interarm 
contrasts. 

3. Discussion 

The existence of leading spiral patterns in a differentially 
rotating disk is a complex issue. Questions such as how 
such patterns can survive the shearing of a disk differential 
rotation come to the fore. How can one infer the presence 
of leading and trailing spiral arms if they are not both 
directly seen? 

We have constructed a simple model using the values 
of NGC 4062 from Table The model is very simple in 
the sense that we do not give a radial dependence of arm 
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Fig. 9. Arm/interarm contrast for the H band image of 
NGC 4062. Open and filled squares are the measures in 
the two arms. The solid line is the mean of the two values 
at each radius. 
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Fig. 10. Same as Figure ||, but for the K' image of NGC 
5248. 



density. The synthetic image has a density equal to unity 
for the main trailing arms, and a density equal to 0.65 for 
the weaker leading patterns (this was chosen to get the 
same ratio \A^\/\A^\). 

The Fourier spectra of these synthetic logarithmic spi- 
rals are shown in Figure [ll]. The synthetic spirals, together 
with the contours for the m=2 component are shown in 
Figure [ij. Note that the leading patterns do not appear 
directly on the contours. Nevertheless, their existence can 
he inferred from the interference patterns. As expected, 
the contours show maxima at the intersection of the two 
(trailing and leading) patterns. 

Therefore, if the synthetic spirals represent incoming 
and outgoing spiral density waves on an axisymmetric 
disk, an interference pattern will be established where stel- 




-40 -20 20 40 -40 -20 20 40 
P P 



Fig. 11. Fourier spectra for synthetic spirals constructed 
with the parameters of NGC 4062 (Table 2). The strong 
harmonics at m=4 and to=6 came because we have no dis- 
persion of the 'density' around the mathematical logarith- 
mic spirals (see discussion in Puerari and Dottori 1992). 
The synthetic spirals are exactly bisymmetrical, and thus 
the amplitude of the odd modes is zero. 

lar density will be larger at locales of constructive inter- 
ference in our near-infrared images. 

Although the phases of the two patterns are differ- 
ent (see Table |^), the maxima are separated by approx- 
imately 90°. To quote Toomre ( |l98lD , "The 90° spacing 
of their successive density maxima [in his models] argues 
eloquently for the presence of trailing and leading waves of 
very similar wavelengths, 

4. Conclusions 

A two dimensional Fourier analysis provides a robust 
method for detecting both trailing and leading spirals in 
galaxies far more distant than M51 or M81. 

We have applied our method to near-infrared images 
of NGC 4062 (optically flocculent) and NGC 5248 (opti- 
cally grand design). In both cases, we have inferred the 
existence of dominant m=2 trailing and secondary m—2 
leading spirals in the Fourier spectra. The consequence of 
two peaks with different winding sense in the m—2 com- 
ponent directly implies spiral arm modulation. 

In each case, the pitch angle of both trailing (T) 
and leading (L) waves is almost the same (P-'"=19° and 
P^=-24° for NGC 4062, and P^=28° and P^=-30° for 
NGC 5248). This was the case also for M81 (Elmegreen et 
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Fig. 12. The synthetic logarithmic spirals, and the result- 
ing contours for the to=2 mode. The galaxy is assumed to 
rotate in a clockwise direction. The dark spirals represent 
trailing arms, while leading arms are indicated by spirals 
in light grey. Note the striking similarity between the con- 
tours in this Figure and those in Figure ^. The existence of 
the leading patterns can only indirectly be inferred from 
the modulation in the spiral arms. 



angles comparable to those of the trailing waves, indi- 
cates that the outgoing waves also have short- wavelengths, 
and this is consistent with the expected group veloc- 
ity of outward-moving, leading waves. This result im- 
plies that the bulge region reflects incoming waves (no 
change of pitch angle), but changes their sense of wind- 
ing. Moreover, the amplification mechanism, which is al- 
ways at co-rotation, must be WASER type II, rather 
than WASER type I. The WASER I mechanism involves 
outward-moving, lon(7-wavelengths, trailing waves formed 
by refraction near the bulge and amplified at co-rotation 
witho ut a change in winding sense, i.e. by superreflection 
(Mark |l976l |1977| ). 

Modal analysis of the structures found here, using also 
the observable velocity dispersions of the stars, could in 
principle determine the disk and halo mass distributions 
and the pattern speeds of the spirals, as was done for M81 



(Lowe et al. 1994). Spiral arm amplitude variations have 



the potential to become a powerful constraint for the study 
of galactic dynamics. 

NGC 4062 represents the very first detection of spiral 
arm modulation in the evolved stellar disk of an optically 
flocculent galaxy. 
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al. 1989). This is highly suggestive that the incoming and 
the outgoing wave-trains have similar wavelengths. 

The amplitude ratio jA'^l/IA-'"! is about 0.5 in both 
NGC 4062 and NGC 5248, wherein the higher amplitude 
is consistently assigned to the trailing mode. The ampli- 
tude ratio is independent of absolute magnitude (-19.44 
for NGC 4062 and -21.19 for NGC 5248) and arm class. 

The arm/interarm contrast increases for both galax- 
ies, but not in a monotonic way. The sinusoidal behaviour 
(seen in the modulation of the intensity) betrays the in- 
terference between incoming and outgoing density waves. 

This study has also demonstrated the efficiency of 
near-infrared images for understanding the mass distribu- 
tions in galaxies which appear quite dis-similar optically, 
but which have much in common when examined in the 
infrared regime. 

Our observations of spiral arm amplitude modulations 



supports the idea originally proposed by Lin (1970) that 
density waves turn around by reflection or refraction inside 
galaxies, presumably in the inner regions where the incom- 



ing, short- wavelengths, trailing waves (Toomre 196£) meet 



the kinematically hot stellar bulge. In the present study, 
the observed symmetry of the leading waves, with pitch 
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